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ryj ] Abstract. 

We have investigated critical phenomena in spin glasses _R :c Yi_ a; Ru2Si2 (R = Dy, Tb, Gd). 
These compounds, where the magnetic moments of rare-earth ions interact by the long-range 
Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction via conduction electrons, has uniaxial 
magnetic anisotropy. The separation of the zero-field-cooled and field-cooled magnetization was 
found only along the c-axis in all compounds, and hence, they are classified into the long-range 
>-^j , Ising spin glass. The magnetic anisotropic energies in these compounds are different from each 

<**] ' other in two orders of magnitude, from 330 K to 1.8 K, however, the critical exponents are 

Q \ similar. It clearly indicates a presence of the universality of the long-range RKKY Ising spin 

O . glasses. 
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1. Introduction 

Spin glass (SG) state is an archetype of ordered states in random magnets where the 
ferromagnetic (positive) and antiferromagnetic (negative) spin-spin interactions are randomly 
distributed. The strong frustration due to the random distribution of the interaction produces 
a complicated magnetic ordered state, where the spins are frozen with random orientations, 
showing striking glassy behaviors such as the long-time relaxation, aging phenomenon, memory 
effect, and chaos effect [H[2]. Although extensive investigations of the SG have been done since 
its discovery, many unresolved problems remains. One of them is whether the replica-symmetry- 
breaking (RSB) predicted by the mean-field theory [3] occurs in real SG materials. 

Two alternative pictures of the SG state are theoretically predicted; the mean-field RSB 
picture and the droplet RS picture. According to theoretical studies of the mean- field model, 
the SG phase transition is a 2nd order phase transition associated with the RSB [?,[3]. The 
RSB SG state has a complicated multi- valley structure of its free-energy [5] and is an essentially 
different state from RS states such as paramagnetic and (anti)ferromagnetic states, where simple 
one- and two-valley structures, respectively. Alternatively, the RS SG state is predicted by the 
droplet theory based on the short-range interaction model [6]. In the droplet picture, the SG 
state has a simple two- valley structure of its free energy as well as in an (anti) ferromagnetic 
state. The RSB and RS SG states can be distinguished by examination of stability of the SG 
phase under finite magnetic field. The RSB SG state is stable in the presence of finite magnetic 
field [7], whereas the droplet RS SG state is unstable in any magnetic field [6]. 



Many numerical and experimental investigations were performed to examine the stability 
of the SG under magnetic field [H O [101 EI]- The experiments of two different model SG 
materials, Fe^Mni-^TiOa (9J [10] and Dy z Yi_a;Ru2Si2 [UJ, are worthy of attention. The 
former Fe^Mni^TiOs and the latter Dy r Yi_ ;! ;Ru2Si2 are model magnets of the Ising SGs 
with the short-range and long-range interactions, respectively. The short-range super-exchange 
interaction via oxigen-ion is dominant in Fe^Mni^TiOs and the long-range Ruderman-Kittel- 
Kasuya-Yoshida (RKKY) interaction via conduction electrons is dominant in Dy x Yi_ :r Ru2Si2. 
The experimental examinations on the stability of the SG under magnetic field were performed in 
both materials to verify critical divergence of a characteristic relaxation time r(T, H) with finite 
SG transition temperature T g (H). In Fe x Mni_ x Ti03, t(T,H) exhibits the critical divergence 
only at zero field [21 [TO] . On the other hand, the critical divergence of r(T, H) in Dy a; Yi_. r Ru2Si2 
was found in both zero and finite fields These results indicate that the RS and RSB SG 
states are realized in the short-range and long-range Ising SGs, respectively. Furthermore, the 
static and dynamic critical exponents of Fe^Mni-^TiOa and Dy :r Yi_ a; Ru2Si2 are quite different 
from each other, i.e., 7 = 4.0, /3 = 0.54, 5 = 8.4, and zv = 11 in Fe^Mni^TiOa [91 [TOl [T2] 
and 7 = 1.1, P = 0.9, 5 = 1.7, and zv = 2.1 in Dy a Yi_. E Ru2Si2 [UJ, indicating that these SGs 
belong to different classes each other and a spacial-range of interaction is relevant to SG phase 
transitions. The critical exponents of Dy r Yi_. E Ru2Si2 are similar to those of the mean-field 
model, i.e., 7mf = 1, /3mf = 1, <5mf = 2, and {zv)mf = 2 [13]. 

In this article, we have reported on static critical phenomena in spin glasses -R x Yi_ x Ru2Si2 
(R = Tb, Gd) in addition to Dy x Yi_ x .Ru2Si2. The parent compounds of these SGs, DyRu2Si2, 
TbRu2Si2, and GdRu2Si2, are isostructural intermetallic compounds, the tetragonal ThCr2Si2 
structure, with paramagnetic-antiferromagnetic transitions of 29 K, 53 K, and 47 K, respectively 
[T4] . In TbRu2Si2 and GdRu2Si2, the long-range RKKY interaction is dominant and spin glass 
(SG) phase also appears by substitution of non-magnetic Y for R, as well as in DyRu2Si2- 
Although magnetic anisotropic energies in these compounds are different in two orders of 
magnitude , from 330 K to 1.8 K, the magnetic anisotropy is the same, an Ising-like, and 
static critical exponents, 7 and 5, are similar to each other. It clearly indicates a presence of 
the universality of the long-range RKKY Ising SGs. 

2. Experimental details 

The single crystalline samples of i? x Yi_ a: Ru2Si2 (R = Dy, Tb, Gd) were grown by the 
Czhochralski method with a tetra-arc furnace. The concentrations of R were determined by 
comparing the saturated magnetizations of the diluted compounds along the magnetic easy c- 
axis with those of the pure compounds i?Ru2Si2- Ac and dc magnetizations were measured with 
the SQUID magnetometer MPMS (Quantum Design) equipped in the Research Center for Low 
Temperature and Material Sciences in Kyoto University. The dc magnetization was measured 
10 minutes after stabilization at a certain temperature and a certain field to avoid the effect 
of very slow dynamics and to obtain the thermodynamic equilibrium magnetizations near the 
SG transition temperatures. The ac measurement of Gd x Yi_ :r Ru2Si2 was performed with an 
ac field of 1 Oe and a frequency of 0.5 Hz. We use thin plate-shaped samples with weights of 
about 7.0 mg and sizes of 5 x 1 x 0.2 mm 3 for the measurements to inhibit the effects of the 
diamagnetic field and the eddy current. 

To examine the static critical behaviors, we extracted the second nonlinear susceptibility 
X2(T) from the dc magnetic field (H) dependences of the dc magnetization M(H,T) and the 
real part of the ac susceptibility x'(H,T). From a dc magnetization measurement, X2(T) can 
be obtained by fitting M(H,T)/ H at each temperature as a function of H 2 in the form of 

M(H,T)/H = Xo (T) + X 2(T)H 2 + ■ ■ ■ . (1) 

X2(T) can also be obtained from an ac susceptibility by fitting x'{H,T) as a function of H 2 in 



the form of 



X '(H,T) = Xo(T) + 3x' 2 (T)H 2 + ---, 



(2) 



because x'(H,T) with small ac field represents a differential magnetization dM/dH. When an 
ac- frequency uj is so low that we can recognized the measurement to be in the dc- limit, X"i(T) 
obtained from the ac measurement is equivalent to X2(T) obtained from the dc measurement. In 
the present work, we obtained X2(T) , s of Dyo.io3Yo.897Ru2Si2 and Tbo.o65Yo.935Ru2Si2 from the 
dc measurements because they show strong frequency dependence of the ac susceptibilities even 
very low frequency of 0.1 Hz. Gdo.o57Yo.943Ru2Si2 hardly shows the frequency dependence and is 
within the dc limit below 1 Hz, and hence, we obtained its X2(T) from the ac measurement with 
u = 0.5 Hz. The field-depedent nonlinear susceptibility Xn\(H,T) = Xo(T) — M(H,T)/H was 
extracted by using the linear susceptibility Xo{T) obtained from above-mentioned procedures. 

3. Results and discussions 

Figure Q] shows T dependences of dc magnetizations of (a) Dyo.io3Yo.897Ru2Si2, (b) 
Tbo.o65Yo.935Ru2Si2 and (c) Gdo.o57Yo.943Ru2Si2 along the c- and a-axes with the applied 
magnetic field of 2 Oe. All compounds exhibit distinct separations of the zero-field-cooled (ZFC) 
and field-cooled (FC) magnetizations along the c-axis around 1.9 K (Dy), 2.7 K (Tb), and 2.0 K 
(Gd). In contrast, no separations of the ZFC and FC magnetizations along the a-axis are found. 
These results imply that only the c-component of magnetic moments of the magnetic Dy 3+ -, 
Tb 3+ -, and Gd 3+ -ions are randomly frozen and the a-component dynamically fluctuates even 
at the lowest temperature of measurements. These component-separated freezing behaviors, 
which was theoretically predicted in the anisotropic SG based on the mean-field model |15j . 
indicate that these compounds should be categorized into the Ising SG with the long-range 
RKKY interaction. 
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Figure 1. T dependences of dc magnetizations of (a) Dyo.io3Yo.897Ru2Si2, (b) 
Tbo.o65Yo.935Ru2Si2, and (c) Gdo.o57Yo.943Ru2Si2. Upper and lower figures are the data with 
magnetic field along the c-axis and a-axis, respectively. Arrows indicate the SG transition 
temperatures. 



Figure [2] shows T dependences of second nonlinear susceptibilities X2(T) , s, corresponding to 
the order-parameter susceptibility of the SG, (a) Dyo.io3Yo.897Ru2Si2, (b) Tbo.o65Yo.935Ru2Si2 
and (c) Gdo.o57Yo.943Ru2Si2 along the c-axis. In all compounds, clear divergent behavior of 
X2(T) are found. From the log-log plots of — X2 vs e (= T/T g — 1) shown in the lower figures 
of Fig. [2j we estimated the critical exponent 7 of —XziX) x e_7 - The best plots showing 
linear relations between log(-X2) and loge with T g = 1.925 K (Dy), 2.740 K (Tb), and 2.008 
K (Gd). The estimated 7's are 1.08, 0.90, and 1.24 in Dyo.i03Yo.897Ru2Si2, Tbo.o65Yo.935Ru2Si2 
and Gdo.o57Yo.943Ru2Si2 respectively. These values of 7 are similar to each other and are 
approximately the value of the mean-field model 7mf = 1- 




Figure 2. T dependences of second nonlinear susceptibilities of (a) Dyo.io3Yo.897Ru2Si2, (b) 
Tb .o65Yo.935Ru2Si2, and (c) Gd .o57Yo.943Ru2Si2. Lower figures are the log-log plots of -X2 vs 
e (= T/T g -1). 

Next, we show H 2 dependences of field-dependent nonlinear susceptibilities at T g 
Xn\(H,T g ), corresponding to the SG order-parameter, in Figs [3] (a) Dyo.io3Yo.897Rrt2Si2, (b) 
Tb .o65Yo.935Ru2Si2, and (c) Gd .o57Yo.943Ru2Si2 along the c-axis. In the figures, Xnl y s- H 2 
are plotted in double logarithmic scales. From these plots, the critical exponent 5 of 
is obtained as 5 = 1.70 (Dy), 1.70 (Tb), and 2.16 (Gd). These values of 5 are similar to each 
other and are approximately the value of the mean-field model 5mf = 2, as well as 7. 

The similarity of the static critical exponents 7 and 5 indicates that the SG phase transitions 
of i? x Yi_ x Ru2Si2 should belong to the same universality class, which is different from those of not 
only Fe x Mni_ x Ti03 [12] but also so-called canonical SGs |16[ I17j. The anisotropic behaviors 
shown in Fig. [1] clearly indicate the Ising nature of i? x Yi_ :r Ru2Si2. We roughly estimated 
magnetic anisotropic energies E^s from the differences between magnetization processes along 
the c- and a-axes; E a = 250, 330, and 1.8 K/i? 3+ -ion in Dyo.io3Yo.897Ru2Si2, Tbo.o65Yo.935Ru2Si2 
and Gdo.o57Yo.943Ru2Si2 respectively. Moderate magnitude of E a (« T g ), which is enough large 
anisotropy to show the Ising nature of SG phase transition according to the theory [15] . is found 
even in Gd-compound with zero orbital angular momentum (L = 0). It should be noted that 
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Figure 3. Log-log plots of H 2 dependences of field-dependent nonlinear susceptibilities of (a) 
Dyo.i03Y .897R.u 2 Si 2 , (b) Tb .o65Y .935Ru 2 Si2, and (c) Gd .o57Yo.943Ru 2 Si2 at each T g . 

the critical exponents in i?. c Yi_; C Ru2Si2 are almost identical in spite of the large difference of 
E a , and are similar to those of the mean-field model. 

Here, we discuss on the validity of our static scaling analyses. The cusp- like anomalies of 
the ZFC magnetization of the Dyo.io3Yo.897Ru2Si2 and Tbo.o65Yo.93sRu2Si2 are rather broad, 
as shown in Figs [T] (a) and (b). It should be a dynamical effect because very slow dynamics 
is found in both compounds [llj. The slow dynamics in the Dy- and Tb-compounds result 
from the strong Ising anisotropy and the consequent high energy barrier of the spin nipping. 
Due to this dynamical effect, it is very difficult to observe the static critical properties very 
close to T g in the Dy- and Tb-compounds. Thus, the regions of the measurements in these 
compounds shown in Figs. [2] and [3] are rather far from T g , as comparing with that in the earlier 
work of the canonical SG AgMn [16]. The analysis of the data in rather high temperature 
and high field may lead wrong critical exponents. On the other hand, the magnetic anisotropy 
of the Gdo.o57Yo.943Ru2Si2 is weak and the effect of the slow dynamics is negligible above T g , 
and hence, we can approach very close T g , as shown in Figs. [2] (c) and [3] (c). Nevertheless, 
we obtained the critical exponents, 7 = 1.24 and 5 = 2.16, which are roughly identical to the 
exponents of the Dy- and Tb-compounds and are much closer to the mean-field exponents than 
to those of the short-range Ising SG or the canonical SG (TB]. Consequently, we concluded that 
-RaYi_ x .Ru2Si2, the long-range RKKY Ising SG, belongs to the mean-field universality class. 

4. Conclusion 

We examined the static critical phenomena of the SG phase transition in i? x Yi_ x .Ru2Si2 (R 
= Dy, Tb, Gd), which are categorized as the long-range RKKY Ising SG, and found almost 
identical critical exponents 7 ~ 1 and 5 ~ 2. These values are similar to the mean-field values 
7mf = 1 and 5mf = 2. It indicates a presence of the universality of the long-range RKKY Ising 
spin glasses, being the mean-field universality class. 
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